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CO ' Abstract. We consider a many particle quantum system, in which each particle 

interacts only with its nearest neighbours. Provided that the energy per particle has 
an upper bound, we show, that the energy distribution of almost every product state 
becomes a Gaussian normal distribution in the limit of infinite number of particles. 
We indicate some possible applications. 
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Physical systems, composed of interacting identical (or similar) sub- 
systems appear in many branches of physics. They are standard in 
condensed matter physics. 

Assuming that each subsystem only interacts with its nearest neigh- 
ed . bours and that the energy per subsystem has an upper limit, which 
must not depend on the number of subsystems n, we show that the dis- 
tribution of energy eigenvalues of almost every product state converges 
to the Gaussian normal distribution in the limit of infinitely many 
^ ' subsystems. To the best of our knowledge, this fundamental quantum 
■ feature has not yet been recognized in the literature [1, 2, 3]. 

Central limit theorems for the distribution of energy eigenvalues 
in quantum gases with Boltzmann statistics [4] as well as for Bose 
and Fermi statistics [5] have been discussed by M. Sh. Goldstein. His 
theorems apply for mixed states, namely classical mixtures of quantum 
states involving classical probabilities. 

We consider here the distribution of energy eigenvalues for a pure 
quantum state. In our case, the distribution is thus of purely quantum 
nature, it solely exists because the state we consider is not an eigenstate 
of the energy operator of our system. Our theorem may be viewed as 
a central limit theorem for mixing quantum systems. 

© 2008 Kluwer Academic Publishers. Printed in the Netherlands. 



hartmann.tex; 7/02/2008; 14:33; p.l 



2 



Some extensions of the central limit theorem to quantum systems, 
with the state not necessarily being mixed, have been proven in the 
past [6, 7, 8, 9, 10, 11]. The version, which appears closest related to 
ours, has been published by Goderis and Vets in 1989 [9]. They consider 
a quantum lattice system and assume that the state and the operator, 
they look at, are invariant under lattice translations. Their proof is 
then based on a set of "cluster conditions", which replace the mixing 
condition of the random variable case. 

The assumptions we use are stricter with respect to the mixing 
behavior, nevertheless, they may still be weakened and generalised. 
On the other hand, we do not assume translational invariance of the 
operators or the state. Instead, we use a quantum analogue of the 
Lyapunov condition for random variables. This generalisation opens 
up a large field of applications. The version of Coderis and Vets only 
applies to products of identical subsystem states, while ours applies 
to almost every product state with the fraction of exceptions being 
negligible. 

Knowing the energy distribution of a product state one can deduce 
estimates on various quantities of interest. If the system is known to 
be in a product state at some initial time, one can calculate it's en- 
ergy distribution. As this distribution is conserved under Schrddinger 
dynamics, one can then make predictions on the dynamics of the state, 
even in the long time limit. These circumstances should prove helpful in 
many problems related to chaotic or non-chaotic behavior in quantum 
systems. 

On the other hand, for a given global state of the total system being 
a function of the total Hamiltonian H, one can calculate occupation 
probabilities of product states. Since only in the product basis traces 
over single subsystems can be performed, this procedure allows to cal- 
culate properties of the reduced density matrix of any subsystem. In 
this way, one could analyse the local properies of stationary global 
states. 

Finally, our theorem may underly the often-used assumption of Gaus- 
sian fluctuations. 

Notation: We first give some notation and definitions we use: We 
consider a chain of quantum systems with next neighbour interactions. 
Let the entire system be described by a Hamiltonian H which is a 
linear, self-adjoint operator on a seperable, complex Hilbert space 
The Hilbert space J4? is a direct product of the Hilbert spaces of the 
subsystems, 



n 




11=1 
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and the Hamiltonian may be written in the form, 

H=J2 n ^ ( 2 ) 

with 

= I®"" 1 ® # M (8) I®" - ^ + ® (8) (3) 

where if M is the proper Hamiltonian of subsystem /i, and i^+i the 
interaction of subsystem /i with subsystem fi + 1. I is the identity 
operator. We chose the boundary condition i^n+i = 0. 

Let E v be the eigenenergies and, using the Dirac notation [12], let 
{\(p)} be an orthonormal basis of Jtff consisting of eigenstates of the 
total system. 

H\<p) = Ey,\<p) with {<p\<p') = S w > , (4) 

where b w < is the Kronecker delta. We denote by \a) the product state 

n 

|a) = n>K>> ( 5 ) 
n=i 

built up from some state |a M ) of each subsystem /z, |a^) S Jif^. 
We furthermore define, 

£ a = (a\H\a) (6) 
cj2 = (a\H 2 \a) - (a\H\a) 2 , (7) 

and introduce the operator 

Z n = (8) 

which is diagonal in the same basis as H. Let denote its eigenvalues, 

Z n \f) = Zp\<p). (9) 

Note that H and therefore E a , a a and z v as well as the basis {\<p}} 
depend on n. 

Since H and thus Z n are self-adjoint, \a) induces a measure on the 
spectrum of Z n respective H. This measure of the quantum mechanical 
distribution of the eigenvalues of Z n in the state |o), is given by the 
usual formula, 

Pa(^e[z 1 ,z 2 ])= K"^)! 2 ' ( 10 ) 

{\lfi).Zl<Z v <Z2} 
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where the sum extends over all states \(p) with eigenvalues in the 
respective interval. 

Theorem: If the operator H and the state \a) satisfy 

a 2 a >nC (11) 

for all n and some C > and if each operator 7i a is bounded, i.e. 

(X|^IX)<C (12) 

for all normalised states \x) G ^ and some constant C, then the 
quantum mechanical distribution of the eigenvalues of Z n in the state 
| a) converges weakly to a Gaussian normal distribution: 

lim F a ( Z(p e [ Zll z 2 \) = r ex P(-_f, 2 /2) dz (13) 

n -*°° J Zl V2-7T 

for all — oo < z\ < zi < oo. 

Proof: Following the proof of the central limit theorem for mixing 
sequences [13] as a guideline, we prove the statement (13) in three steps: 
First, we show that the characteristic function of H does not change if 
a few of the TL^ are neglected. Second, we prove, that the characteristic 
function of the remainder of H factorises. In the last step, we then show 
that the condition for Lyapunov's version of the central limit theorem 
is fulfilled for the remainder of H. The proof is then completed by 
the standard proof of Lyapunov's central limit theorem, which can be 
found in several textbooks [14]. 

Define the operators = Ti.^ — (a|?^|a) and split the sum 

1 n 

Z n = - £ X, (14) 
° a ,1=1 

into alternate blocks of length k — 1 (large blocks) and of length 1 (small 
blocks). The large blocks are given by 

k-l 



£j = J2 X {j~i)-k+i for 3 = !> • • • » \ n / k \ and ( 15 ) 
l=i 

Q 

f[n/fc]+i = zZ X lWk]-k+i with q = n — k[n/k], (16) 

where [x] means the integer part of x and the small blocks are the 
Xj.k with j = 1, . . . , [n/fc]. Sum up all large blocks and all small blocks 
separately, 

[n/k] + l [n/k] 
Z 'n = - E & and Z n = ~zZ X 3*, (17) 

°a j=1 a a j=l 
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so that Z n = Z' n + Z'n- The integer block length k is chosen to depend 
on n (k = k(n)) such that 



n k 
hm —7; = and lim — = 0, 

n— >oo k n—>oo ji 



(18) 



with k 



n 



3/4 



being a possible realisation. 
: fun 

(a\e~ 



Consider the characteristic function 

-irZ n 1 



with real r. 

First let us show that for all real r and n — > 00: 



(a\e- irZ "\a) (a|e 



Using the operator identity [15] 

Jo 

the triangle- and the Schwarz-inequality, one gets 



(19) 



(20) 



(21) 



(a\e- irZn - e~ irZ " 



a)\ < I dsJ (a\e isZ " (Z%) 2 e- isZ "\ 



< r 



\ 





~n~ 


2 \ {2C'Y 


\ n 


k. 


J c 



(22) 



which, indeed, converges to zero for n — > 00. 

Next, we show that the characteristic function of Z' n factorises. 



(a\e 



-irZL 



[n/k]+l 

n ( a \ e ~^\ 



(23) 



To this end, we first note two important properties that arise due to 
the next neighbour interaction and the product property of the state 
I a): For |/i — u\ > 1 and any two integers k and I, we have 

[n^Hu] = (24) 
(a\ {H^f {H u ) [ \a) = {a\{H^) k \a){a\(H u ) l \a). (25) 

Therefore, for all and any two integers k and I, 

[Si, €3] = (26) 
(a\&) k (& l \a) = (a|&)*|a>(a|&)'|a>, (27) 
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and equation (23) follows as a direct consequence. 

Finally we prove that the £j fulfill the Lyapunov condition: 



1 



[n/k]+l 

E 



1™J^ E (a\\^\ 2+m \a) = (28) 



for some m > 0. Note that due to equation (20), (a| {Z' n ) 2 \a) — > a 2 a as 
n —> oo and therefore equation (28) is, indeed, the Lyapunov condition 
for the £j. We verify the condition for m = 2. To this end, consider 

fe-i 

(a|^|a) = ^2 (a\X(j_i) k+ pX(j_ 1 j k+v X(j_ 1 j k+p X(j_i} k+T \a). 

ll,V,p,T = l 

(29) 

Since (a\X^\a) = and because of equations (26) and (27), only those 
terms are nonzero, for which all the X^ are identical or neighbours or 
where two pairs of identical or neighbouring Xn appear. For example 
{a\X^X^ + iX^ +2 X^i\a) / while (a\X ^X M+3 X M _i|a) = or 
(a\X^ 1 X il+ iX u ^iX v \a) / while {a\X fl X u X^X u+2 \a) = 0. Using this 
fact and the conditions (11) and (12) one realises that 

[n/fc]+l 

~, E «l«> < 



< 



+ 1) «*-!)'+ 3-5.7 ^- (»-!)) (SOT (3Q) 



where the rhs vanishes in the limit n — > oo. Note that ( a l£l/fc] + il a ) 
contains less terms than (a|£||a) for j < [n/k] + 1 and is therefore 
bounded by the same expression. 

With the arguments above showing that the characteristic function 
of Z' n factorises and that the £j obey the Lyapunov condition, it is 
straight forward to prove, following the standard steps [14], that 

( r<2 \ 

lim (o| exp (—irZ ) \a) = exp — — I (31) 

n^oo v ' \ 2 I 

and, using equation (20) one concludes that 

( r2 \ 

lim (a exp (—irZ n ) la) = exp — — 1 . (32) 

n^oo \ 2 1 

Here both limits are pointwise for all real r. 

Finally, the continuity theorem [14] states that the pointwise con- 
vergence of the characteristic functions, established above, implies the 
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weak convergence of the distributions. The density of the limit dis- 
tribution is thus given by the Fourier transform of the characteristic 
function in equation (32), which proves our theorem. 

Note that all the steps not explicitely carried out here only use prop- 
erties of Lebesgue integration and no further properties of probability 
distributions. We therefore do not run into difficulties related to so 
called "no hidden variable theorems" [16, 17]. 

For applications in physics, where n is very large but finite, the 
density of the limit distribution can be written as a function of the 
energy E of the system, 



Pa{E) = -= — exp 
V27T a a 



V a J 



(33) 



so that P a (E G [E 1 ,E 2 ]) = pa(E)dE. 

Discussion and Generalisations: Let us first analyse the condi- 
tions (11) and (12) in more detail. 

Rewriting (11) in terms of the operators and using equation (25) 
we get 

n -. 

Y,(a\-(xl + Xl +1 ) + X^X^^ X^X^a) >nC. (34) 
,u=l 

Thus, every term in the sum in (34) being larger than C is sufficient 
for (11) to be satisfied. 

Condition (12) physically states, that the excitation energy must 
not be concentrated in only a small part of the subsystems. For very 
large systems, where our theorem applies, this is only a minor restric- 
tion since the fraction of states that do not fulfill condition (12) is 
vanishingly small. 

Several conditions we have used to derive our theorem may be 
relaxed and substituted by weaker assumptions. 

First of all, the theorem is not only valid for a linear chain but also 
for two and three dimensional lattices. 

It is also straight forward to proof the same theorem for periodic 
boundary conditions, I n ,n+i = Inf- 
ill addition, it is obviously not necessary that the subsystems only 
interact with their nearest neighbours. The theorem holds as long as 
the number of interaction partners (the connectivity) of each particle 
is limited. 

Furthermore, the observable one considers need not be the Hamilto- 
nian. Any other observable shows the same feature as long as conditions 
(11), (12), (24) and (25) are met. 



hartmann.tex; 7/02/2008; 14:33; p. 7 



8 



Let us stress here, that neither the operator (Hamiltonian) nor the 
product state need to be invariant with respect to lattice translations. 

Finally, conditions (11) and (12) may be relaxed, since the theorem 
still holds, whenever Lyapunov's condition, or even only Lindeberg's 
condition [14], is fulfilled. We have chosen here stricter but simpler 
conditions to make it easier to check the applicability of our theorem. 

Examples: We mention two examples, where our theorem applies. 
First, we consider an Ising spin chain of the type [18] 

ff = -B5>?- | Safari. (35) 

i i 

Here, sf and sf are the Pauli matrices, 

11) md s * = U-i)- (36 > 

B the difference between local energy levels and J the coupling strength. 
The energy of each spin is at least —B and at most B so that condition 
(12) is satisfied. The squared width a 2 a reads 

*a = «X' (37) 

where n is the number of spins, and condition (11) is also met. 
As a second example, we consider a harmonic chain. 

# = -£f| + ^ 2 te+i-<fc) 2 (38) 

i 

where qi and pi are the position and momentum of particle number 
i. All particles have mass m and the coupling has frequency u. Since 
the energy of a harmonic oscillator is not bounded, our theorem only 
applies to states where the energy per oscillator does not exceed a 
certain bound, which on the other hand may be chosen arbitrarily large. 
The squared width for the harmonic chain is 

°l = nu? (ni + ^j + (39) 

where rtj is the occupation number of oscillator number i. Since rtj > 0, 
condition (11) is satisfied. 

Applications: We finally discuss two areas of possible applications 
of our result. 

a: One may consider the product state \a) as an initial state and 
make predictions about its dynamics. Using equation (32), one can 
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calculate the fidelity of state |o): 



-iHt 



2 




(a\e 



a) 



(40) 



Furthermore, one can give an upper bound to the transition probability 
to another product state \b) ({a\b) = 0) for all times t; 



where we have assumed a a <C E a and <C E\, and the ground state 
energy has been chosen to be zero (Eq = 0). 

b: One can calculate diagonal elements of the reduced density matrix 
of a selected subsystem in the basis |o M ), provided the total system is 
in a stationary state, that is, its density matrix ptotai is a function of 
H [19]. For this applications, it is most interesting to take \a) to be 
an eigenstate of the Hamilton operator without the nearest neighbour 
interactions (Ho\a) = E a \a), Hq = Y^uHp). Here, very interesting con- 
clusions can be drawn on the minimal spatial extension of temperature, 
which will be presented elsewhere [20]. 

In summary, we have considered a large quantum system composed 
of subsystems, where each subsytem only interacts with a limited num- 
ber of neighbours. We have shown that for almost every product state, 
the distribution of the total energy converges to a Gaussian normal 
distribution in the limit of infinitely many subsystems (13). This is the 
main result of this paper. The assumptions we have made are quantum 
mechanical analogues to the conditions for Lyapunov's central limit 
theorem for mixing random variables. We did not dwell on the most 
general assumptions needed for our theorem to make the verification 
of our conditions in physical applications straightforward. Nevertheless 
we have discussed possible generalisations as well as some preliminary 
applications. 

We thank J. Gemmer, M. Michel, H. Schmidt, M. Stollsteimer, F. 
Tonner, M. Henrich and C. Kostoglou for fruitful discussions. 

M.H. in particular wants to thank Prof. Detlef Diirr for many help- 
full comments. 



1. Licb, E. and Mattis, D.: 1966, 'Mathematical Physics in One Dimension', 
Academic Press, New York 

2. Thirring W.: 2002, 'Quantum Mathematical Physics' Springer , Berlin, 2nd ed. 




(41) 



References 



hartmann.tex; 7/02/2008; 14:33; p. 9 



10 



3. Lieb, E.: 2001, 'The Stability of Matter: From Atoms to Stars', Springer , Berlin 

4. Goldstein M. Sh.: 1985 J. Stat. Phys. 40, pp. 329 

5. Goldstein M. Sh.: 1985 Theoret. and Math. Phys. 1, pp. 412 

6. Cushen, CD. and Hudson, R.L.: 1971 J. Appl. Probability 8, pp. 454 

7. Accardi, L. and Bach, A.: 1985, Z. Wahrsch. verw. Geb. 68, pp. 393 

8. Goderis, D., Verbeure, A. and Vets, P.: 1989, Prob. Th. Rel. Fields 82, pp. 527 

9. Goderis, D. and Vets, P.: 1989, Comm. Math. Phys. 122, pp. 249 

10. Kuperberg G.: 2002, math-ph/ 0202035 

11. Michoel, T. and Nachtergaele, B.: 2003, math-ph/0310027. 

12. Sakurai, J. J.: 1994, 'Modern Quantum Mechanics', Addison- Wesley, Reading, 
Massachusetts 

13. Ibargimov, LA. and Linnik, Y.V.: 1971, 'Independent and Stationary Sequences 
of Random Variables', Wolters-Noordhoff, Groningen/Netherlands 

14. Billingsley, P.: 1995, 'Probability and Measure', John Wiley & Sons, New York, 
3rd ed. 

15. Fick, E. and Sauermann, G.: 1990, 'The Quantum Statistics of Dynamic 
Processes', Springer , Berlin 

16. Diirr, D., Goldstein, S., Zanghi, N.: 2003, J. Stat. Phys. 67 pp. 843 

17. Diirr, D., Goldstein, S., Zanghi, N.: 2004, J. Stat. Phys. to be published, 
quant-ph/0308038 

18. Mahler, G. and Weberrufi, V.: 2001, 'Quantum Networks', Springer , Berlin 

19. Jordan, A.N. and Biittiker, M.: 2003, cond-mat/0311647 

20. Hartmann, M., Mahler, G. and Hess, O.: quant-ph/0312214 



hartmann.tex; 7/02/2008; 14:33; p. 10 



